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Abstract: Gravitational lensing can magnify a distant source, revealing structural detail 
which is normally unresolvable. Recovering this detail through an inversion of the influ- 
ence of gravitational lensing, however, requires optimisation of not only lens parameters, 
but also of the surface brightness distribution of the source. This paper outlines a new 
approach to this inversion, utilising genetic algorithms to reconstruct the source profile. 
In this initial study, the effects of image degradation due to instrumental and atmo- 
spheric effects are neglected and it is assumed that the lens model is accurately known, 
but the genetic algorithm approach can be incorporated into more general optimisation 
techniques, allowing the optimisation of both the parameters for a lensing model and the 
surface brightness of the source. 

Keywords: Gravitational lensing — methods: numerical 
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Gravitational lensing is an important astrophysical tool, 
mapping the distribution of matter on many scales and 
revealing typically unresol ved detail of dista nt sources 
through magnification [see iKochanek et alJ (|i]04) for 
a review of lensing physics]. Early studies concen- 
trated upon the properties of multiply imaged quasars 
to determine the underlying mass distribution within 
the lensing galaxy. However, the point-like nature of 
quasars (which results in point-like images) present 
only a limited amount of constraints on the lensing 
mass distribution, with a number of degenerat e solu- 
tion s able t o explain the observe d configuration iKent fc 
[1^8 ; Schne ider etHIIllQSSl: [KQchanck . 

If the source in a gravitational lensing is extended, 
the resulting image is also extended and each resolu- 
tion element effectively provides a constraint on any 
modelling. Any modelling, however, becomes more 
complex than the simple case of point-like lensing which 
asks the question "what distribution of mass in the 
lensing can account for the observed image locations 
and brightnesses?". With an extended source, the 
question has to be rephrased as "what distribution of 
mass in the lensing galaxy and distribution of bright- 
ness in the source can account for the observed image 
configuration?". To answer this question, more novel 
approaches to gravitational lensing modelling have been 
und ertaken (Kochanck. Blandford. Lawrence. & Narava3| 
1%9; Kochanck & Naravan 1992; EUithorDC, Kochanck, & 



problem, which is what occurs in practice, is much 
more difficult to solve. Furthermore, inverse problems 
are fraught with the question of solution uniqueness. 
To this end, it is advisable to tackle a particular prob- 
lem with a range of inversion techniques and compare 
the various outcomes; if all approaches converge to the 
same result, some faith can be given to the overall solu- 
tion. Currently, the repertoire of gravitational lens in- 
version techniques is relatively small and so this paper 
presents an alternative approach to gravitational lens 
modelling, utilising genetic algorithms to reconstruct 
the source. This initial investigation of this approach, 
th e effects of image degradation through instrumental 
Falcal id atmospheric effects (i.e. seeing) are neglected, al- 



though these can be implemented in a straightforward 
fashion. The approach is described in detail in Sec- 
tion |21 considering a perfect, noiseless image. Section 
13 discusses the influence of the various parameters in- 
fluencing the inversion technique, also considering re- 
construction of noisy images. Section |1] considers the 
more general problem of the optimization of both the 
source surface brightness distribution and the param- 
eters governing the mass distribution in the lensing 
galaxy. This paper closes with the conclusions which 
are presented in Section |3 



200a : lwucknitzl2004l:IWavth et alJ2005l : lDve fc WarrenI 
2003) • These techniques typically determine not only 
the mass distribution in the lensing galaxy, but also 
the surface brightness distribution in the source. 

Given the form of the lens and the source, it is rela- 
tively straightforward to compute the resultant image 
by using a simple ray tracing method. The inverse 



2 Approach 

J Hcwtt Qenetic Algorithms 

Dye ^ 

Genetic algorithms are an approach to problems of op- 
timisation that take their inspiration from evolution- 
ary biology [for a popular review of genetic algorithms 
and o ther aspects of biological computing, see iLevvl 
1I1993I ) ]. The basic approa ch, presented in some de- 
tail in ICharbonneaiJ ijlQQflTl . mimics the evolutionary 
struggle of life, with different individuals having differ- 
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Figure 1: The left-hand panel presents the artifi- 
cial source utilised in this study, while the right- 
hand panel presents the image of this source as 
seen through the model gravitational lens (see sec- 
tion Note that these images as on different 
scales, with the source panel being half the width 
of the image panel. The source is defined on a 
grid of 32x32 pixels, whereas the resultant image 
is 64x64 pixels. For observed gravitational lens 
images, the image would be ~ 3 arcseconds on a 
side, corresponding to a pixel scale of 0.05 arc- 
seconds per pixel. The corresponding source plane 
region possesses an image scale of 0.025 arcseconds 
per pixel and a side length of 0.8 arcseconds. 



ent probabilities of passing on their genes to the next 
generation, with the probabilities dependent on the 
environment and the physical characteristics of the or- 
ganism (which in turn depend on the genes). Humans 
have long used a basic knowledge of heredity to breed 
for desirable traits in animals and crops, even before 
Char les Darwin pr oposed his theory of natural selec- 
tion iParwin 11859 !). Genetic algorithms take this idea 
and apply it by breeding better solutions to the prob- 
lem at hand. 

In terms of the algorithmic approach, several fea- 
tures are required; 

• Encoding: Each potential solution to a prob- 
lem is encoded into a genome. This is often rep- 
resented as a series of digits, but can be a simple 
bit string. 

• Expression: This decodes the information in 
the genome into a phenotype. For many appli- 
cations, this decodes the genome into a series of 
real numbers that are used as parameters for a 
particular model. 

• Fitness evaluation: This compares the phe- 
notype of the genome to the problem, assigning 
a quantitative measure of the goodness of fit of 
the potential solution (e.g. a standard mea- 
surement). 

With these, an initial population of genomes, each rep- 
resenting a potential solution to the problem at hand, 
can be generated. Typically this involves assigning 
each genome with a random sequence of digits or bits. 
In evolving this population to the next generation, sev- 
eral steps are involved: 



a) Ranking & Selection Pressure: The goal of a 
genetic algorithm is to produce subsequent generations 
of solutions with greater fitness by ensuring that the 
fittest member of a current population to pass their ge- 
netic material onto to the next generation. As noted by 
Charbonncau (1995), as evolution proceeds, the aver- 
age and maximum fitness of a population continually 
increases. The spread in fitness, however, decreases, 
with the overall population becoming homogeneous. 
With such uniform fitness in a population, simple se- 
lection on fitness alone eff'ectively samples randomly 
from the population and evolution stalls. To circum- 
vent this, selection must be made relative to the cur- 
rent population. To this end, the population is ranked 
in terms of its fitness, with the least fit being assigned 
a value of 0, whereas the most fit possesses a value of 
1. Members are selected from this current generation 
with a probability dependent upon their ranking, such 
that 

p oc (ranking)^ (1) 

where /3 is called the selection pressure. If /3 = 0, 
the probability for selection is uniform throughout the 
population, while larger values of /3 preferentially se- 
lect only the fittest members in the population. 

b) Elitism: The fittest member of the current gen- 
eration is cloned and represents the first member of 
the next generation. This ensures that the maximum 
fitness of subsequent generations can never fall. 

c) Breeding: Further members of the next generation 
are produced by breeding the members of the current 
generation, with the genetic information of the current 
population used to determine the genomes of the next. 
The probability that an individual will breed is based 
upon its ranking and the selection pressure as outlined 
previously. Two breeding strategies are adopted, asex- 
ual and sexual reproduction. With asexual reproduc- 
tion, a selected member of the current generation is 
cloned, preserving the genetic information, to provide 
a new member of the next generation. In sexual repro- 
duction, two members of the current generation are 
selected and a new individual is formed with the com- 
bination of their genetic material; a random portion of 
the genome of one parent is "copied and pasted" over 
the corresponding part of the other parent's genome 
to produce the resulting offspring genome; this allows 
the genetic material of successful organisms to mix. 
Whether the creation of an offspring is due to sexual re- 
production is determined randomly, with a pre-defined 
probability referred to as the crossover rate. 

d) Mutation: A population which reproduces purely 
asexually rapidly becomes dominated by a single genome 
and evolution grinds to a halt. Random mutations in 
the genetic sequence can drive evolution beyond this 
point, increasing diversity in a population ^ The prob- 
ability that a particular digit in the genetic sequence 
is mutated (straight after breeding) is determined by 
the mutation rate. 

The breeding and mutation are continued until a 
new generation is formed, and the entire initial gener- 

^It should be noted that the neo-Darwinistic view that 
it is solely DNA that evolves has been questioned, with the 
implication that evolution is actu ally a complex interplay of 
the genotype and phenotype (see lCohen &: StewartlEoOOfl . 
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Figure 2: An example of the evolutionary se- 
quence obtained with the genetic algorithm recon- 
struction of an idealised gravitational lens system. 
The left-hand panel displays the source plane sur- 
face brightness distribution at 10, 100, 500, 1000 
and 5000 generations, while the right-hand panel 
presents the resulting image plane distribution. It 
is clear, when comparing to Figure ^ that the fit- 
ness of the solution is increasing with each gen- 
eration. Note that the vertical noise at the edge 
of the source plane correspond to regions which 
are not mapped into the image plane and so do 
not contribute to the overall fitness of the image 
reconstruction. 



ation is culled. Steps a) through d) are then repeated, 
with subsequent generations exhibiting fitter solutions 
to the problem. The evolution is terminated when an 
ap propriate fitness crite rion i s satisfi e d. 

ICharbonneeiul \iS9$) and lHakalal \199^ presented 
some of the earliest applications of genetic algorithms 
in astronomy, with the former detailing a freeware 
algorithm (PIKAIA) ^. These early studies focused 
upon t he fitting of light cu rves and galactic rotation 
curves ()CharbonneaxJll99Rr) . and constraining accre- 
tion stream mapping in eclipsing polars (Hakala 1995), 
but since then the apphcation of genetic algorithms 
has been used in the design of filters and filter syst ems 
(Off er fc Bland-Hawthornll998HBailer-Jonesl2004r). mod - 
elling the structure of the ealaxv (|Larsen fc H^mphr^y^ 

2^23) , the signature of g amma ray burst jJPgr^ggig^Zw art fc Totanil 

2001") , solar oscillations (iFletcher. Chaplin, fc Elswor tjj 

200^ and even the scheduling of telescopes (l|G6mez de Castro fc Yane3 

2.2 Gravitational Lens Inversion 

2.2.1 Encoding & Phenotyping 

For the purposes of this project, each individual genome 
was taken to be a string of 1024 characters. The ex- 
pressed genome (the phenotype) represented 32x32 
pixels, where each pixel could take a value between 
and Fmax- This pixel array represents the surface 
brightness distribution of the source. Such an encod- 
ing ensures the surface brightness is subject to a pos- 
itivity constraint [this is not the case of some other 
inversion approaches e.g. IWarren fc Dvel ([2003^ ]. In 
the coming simulations, Fmax = 255. 

2.2.2 Gravitational Lens Model 

The scaled lensing equations relating a position {x, y) 
in the image plane to the corresponding point {xs,ys) 
in the source plane [see IWambseanssI ()l99?^ l] can be 
written in terms of a potential (f>{x,y) (related to the 
mass distribution of the lens) such that 



— _ ££1 

Vs — y ay \ {=',y) 
The gravitational lens potential was chosen to be 
the three parameter pseudo-isothermal elliptic poten- 
tial (Kochanek. Blandford. Lawrence, fc Naravan.l989ll 
of the form 

0(a;, y) = b^rl + (1 - e)x^ + (1 + e)y2 (3) 

where Vc is the core radius of the potential, e is the el- 
lipticity and b is an overall normalisation factor (typ- 
ically linked to the velocity dispersion of the lensing 
galaxy). While the determination of these parameters 
is the typical goal of many analyses of gravitational 
lenses, in this intial examination of genetic algorithms, 
it is assumed that the model is fixed and its parame- 
ters are known. In other words, the goal is to find the 
source profile, given the observed image and the form 

^Jittp : //www . hao . ucar . edu/publlc/research/sl/pikaia| 
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Figure 3: The influence of the mutation rate on 
fitness as a function of generation. The box in this 
panel (and in subsequent figures) denotes the dif- 
ferent values of parameter adopted. Note, this fig- 
ure presents fitness~^ as defined in Equation^ A 
reasonable rate of mutation is needed to stimulate 
progress in the population by providing random 
variations for selection to act upon. However, if 
the mutation rate is too high, any "good" genomes 
will be severely affected by too many mutations. 
The black curve is for the adopted mutation rate 
of IQ-^. 



of the lens. The adopted values for the lens param- 
eters were b — 0.5, e = 0.25 and Vc — 0.1. The full 
optimisation problem of the source profile and model 
parameters will be discussed in Section 0] 

2.2.3 Fitness Determination 

To test the effectiveness of the genetic algorithmic ap- 
proach to gravitational lens inversion, an example solu- 
tion was defined. This consisted of two offset Gaussian 
profiles of differing heights and is displayed graphically 
in the left-hand panel of Figure while the image of 
this source as seen through the lensing potential out- 
lined previously appears in the right-hand panel. The 
image of the source is realised upon a 64x64 grid. 

In determining the fitness of a particular genome, 
its expressed phenotype is used to represent a potential 
source. This is mapped, via the lens model, to produce 
the resulting image configuration. This is then com- 
pared to the ideal image (Figure and the assigned 
fitness was chosen to be the reciprocal of the sum of 
the squared differences between the reconstructed im- 
age and the observed image: 

fitness = ^64^ .._ p (4) 

where rriij are the pixel brightnesses of the model gen- 
erated from the genome and pij are the pixel bright- 
nesses of the "observed" image. Clearly, in the absence 
of noise, the sum of the squared differences for a per- 
fect genome is zero and hence the fitness is infinite. Of 
course, for a real gravitational lens system, this ideal 



Figure 4: As for Figure |3 but detailing the influ- 
ence of cross-over rate on fitness as a function of 
generation. 



image will be replaced with a noisy observed image. 
Simulations of the genetic algorithm considering the 
influence of noise, will be discussed in Section f3. 21 



3 Results 

Examples of the source reconstructions achieved with 
the genetic algorithm are presented in Figure|21 These 
were produced by optimising the source profile, assum- 
ing that the correct lens model was known. From top 
to bottom, each pair of panels presents a snapshot of 
the evolution, showing the fittess member of a gen- 
eration. The left-hand panel graphically presents the 
surface brightness of the source (the genome) , whereas 
the right-hand panel presents the resultant image con- 
figuration (the phenotype). Clearly, as the population 
is evolved to older generations, the accuracy of the 
solution increases (compare the source and image con- 
figuration at Generation 5000 to Figure 0. Note that 
the 'noisy' pixels along the vertical sides of the source 
reconstruction in Figure|5]correspond to regions which 
are not mapped into the image region and so does not 
contribute to the overall fitness of the solution. 

3.1 Evolutionary Parameters 

There are several parameters that can affect the per- 
formance of a genetic algorithm. It is important to 
find a set of parameters which is good at finding a 
solution quickly, since it can be a computationally ex- 
pensive task to run for many generations. Also it may 
be necessary to evolve the population many times, in 
which case it is very important to minimise the time 
required to find a satisfactory solution. There are sev- 
eral factors which have the potential to infiuence this 
performance, and the four which are investigated here 
are the mutation rate, the crossover rate, the size of 
the population and the selection pressure, which is the 
relationship between the fitness score and the proba- 
bility of being selected for breeding. 
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Figure 5: As for Figure 13 but detailing the influ- 
ence of selection pressure on the fitness as a func- 
tion of generation. 
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Figure 6: As for Figure 13 but detailing the influ- 
ence of population size on the fitness as a function 
of generation. 



It is known that there are interactions between 
these parameters such that, for example, the answer 
to the question "what mutation rate should I use?" de- 
pend s on the values of the other parameters ijCharbonneaxJ 
llQQST) . It is also highly dependent on the nature of 
the problem at hand. Therefore, the best that can 
be hoped for is some general picture of what order of 
magnitude these parameters should be. 

To test the dependency of the genetic alorithmic 
approach on the adopted parameters, an initial recon- 
struction was undertaken with the following parame- 
ters; the crossover rate was set to 0.9, the population 
size was 50, the mutation rate was set to 10""^ and 
the selection pressure was set to 10 (these were good 
values found after much painstaking trial and error). 
In the following sections, one of these parameters was 
varied, while the others were kept fixed, allowing the 
influence of the varied parameter to be determined. 



3.1.1 Mutation 

Figure |3 presents the evolution of the "Best Error" 
(defined as fitness'^) of the fittest member of a pop- 
ulation as a function of generation for several differ- 
ing mutation rates. Clearly, if the mutation rate is 
low (10"") then the genome evolves slowly, but shows 
steady improvement. Increasing the rate of mutation 
increases the rate of improvement in the genome. How- 
ever, this cannot continue indefinitely; as can be seen, 
a large mutation rate (lO"'^) results in a rapid increase 
in improvement initially, after a short time the evolu- 
tion stagnates. This is because for a mutation rate of 
0.01 and a genome length of 1024, on average there are 
about 10 mutations per generation on each individual, 
and this is enough to outweigh any improvements that 
have evolved through selection. Hence, a little muta- 
tion is a good thing, but too much mutation is not. A 
general rule of thumb is that a good mutation rate is 
one which will give about 1 mutation over the whole 
genome. 

3.1.2 Cross-over Rate 

Figure 2] presents the results of varying the cross-over 
rate on the rate of improvement of the genome. For 
the three trial values of 0.1, 0.5 and 0.9, there is very 
little difference on the rate of improvement over time. 
Clearly, the algorithm is generally insensitive to the 
adopted value of the cross-over rate. 

3.1.3 Selection Pressure 

Figure |3 presents the inffuence of the selection pres- 
sure, /3, on the rate of improvement of the fittest genome 
Three /3 values of 1, 5 and 10 were trialed; remember 
that /3 = ensures a uniform selection probability for 
breeding from a population (i.e. no selection pressure 
at all) , whereas larger values preferentially selected the 
fittest members for breeding. For the adopted values, 
there was a definite advantage in using a (3 value big- 
ger than 1, but only a slight difference in performance 
between /3 = 5 and P = 10. 

3.1.4 Population Size 

The number of individuals in the population can also 
be varied, and the inffuence of changing population 
size is presented in Figure |S| Clearly larger popu- 
lations have a larger spread in genetic variation and 
it is seen that the larger populations do evolve more 
rapidly. From a computational point of view, however, 
smaller populations result in a significant speed advan- 
tage, with less calculations required per generation. 

In fact, as the time taken for each generation is 
roughly proportional to the population size, the larger 
genetic variability seen in larger population can be out- 
weighing the time required to calculate the fitness of a 
population. 

3.2 Source Reconstruction & Noise 

In real life, astronomical images are always contami- 
nated with some random noise, arising from sources 
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Figure 7: The results of the genetic algorith reconstruction using the noisy image after 5000 generations. 



such as thermal fluctuations in electronics, the efi'ects 
of the atmosphere, and even photon counting noise if 
the source is faint (a feature that is shared by a num- 
ber of extended gravitational lenses). Therefore it is 
worthwhile to test how the genetic algorithm performs 
when the observed image is contaminated with noise. 

Including errors into the observed image also sets 
a meaningful criterion to decide whether a particular 
reconstruction is good enough. This criterion will be 
met if the reconstructed image matches the observed 
image to within the statistical error level set by the 
noise. For this test, a normally distributed random 
variable with a = 5 was added to each pixel in the 
observed image, creating the noisy image shown in the 
lower left-hand panel of Figure 

When the genetic algorithm was run, the error vs. 
time plot (Figure (HJ flattened out at a much higher 
value. This was not surprising, because in this case 
it is not possible to have a source that reproduces the 
observed image exactly, with its random fluctuations 
between neighbouring pixels. The additional panels in 
Figure 13 presents the reconstructed image and source 
profile, as well as the difference between these and the 
true source and corresponding image; these are consis- 
tent with the input noise characteristics. 

This process can be thought of as a curve-fitting 
problem in 2 dimensions, with 64^ data points (the 
observed image), each with an error bar of 5 units. 
The aim is to fit this data using a model that has 
32^ — A'^ free parameters (each pixel of the source; note 
~ 80 corresponds to those pixels in the source grid 
which are not lensed into the flnal image, and hence 
are not true free parameters; see Figure |5J. The chi- 
squared statistic for this flt (with degrees of freedom 
equal to the number of constraints minus the number 



of free paramaters in the model, i.e. Q4^ — 32^ + 80) is 
given by 



64 64 

1=1 i=l 



■Pij) 



■ fitness 



(5) 



Hence, a reconstruction is statistically good (within 
la) if the sum of the squared differences between the 
observed image and the image of the reconstructed 
source is within x {v i: V2i/), where u is the num- 
ber of degrees of freedom. iPress et al.lll993) Hence, a 
good fit corresponds to a value in the range of 78800 ± 
1985. The genetic algorithm was able to reduce the 
error to ~ 82000, and so recovered an acceptable flt to 
the noisy image. 

4 Full Optimisation 

4.1 Lens Parameters in the Genome 

In general, the goal of gravitational lens reconstruc- 
tion is to determine both the surface brightness of the 
source and the parameters describing the mass distri- 
bution in the deflecting galaxy. How can the genetic al- 
gorithmic approach be generalised to tackle this prob- 
lem? One idea that was tried was incorporating the 
lens parameters as part of the potential solutions to 
the problem, by encoding them in the genomes. Then, 
the evolution would hopefully select individuals whose 
lens parameters were close to the right values, and 
then proceed to optimise the source pixels. The re- 
sult of this approach was not very successful. What 
actually occurred was that early on, a particular value 
of parameters was locked on and became dominant in 



www. publish . csiro. au /journals /pasa 



7 




600 1000 1600 2000 2600 3000 3600 

Generation 



Figure 8: Error of the fittest solution vs generation 
for the noisy image. Notice that it does not tend 
to zero, as a perfect reconstruction of this image 
is impossible. 
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Figure 9: The best error after 3500 generations 
as a function of the value of b. Clearly, a good 
reconstruction of the image is only possible if the 
value of b is close to the true value of 0.5. 



the population. Then, the sources were optimised for 
those (incorrect) values, and there was little hope in 
ever getting to the right values. This was because any 
change to the parameter values would require a huge 
chance jump in the sources in order to gain a higher 
fitness than what had already evolved. 

4.2 Direct Search Method 

While incorporating the lens parameters in the genome 
was found not to be successful, an approach using a 
started direct search, independent of the genetic algo- 
rithm, was found to be successful. For the purposes of 
this study, simple one-dimensional searches were em- 
ployed, although the technique can be easily gener- 
alised into higher dimensional searches. 

An example of such a one-dimensional search is 
presented in Figure El in which rc and e were held at 
their optimum value, and a suite of reconstructions 



were undertaken with differing values of b. This fig- 
ure presents the residual of the fittest solution of the 
population after 3500 generations. This clearly reveals 
that a good reconstruction is only possible if b is very 
close to the correct value of 0.5. 

The major disadvantage of this approach is its com- 
putational inefficiency. Each point in Figure|5|required 
the genetic algorithm to be run over 3500 generations, 
which takes about 10 minutes on a modern desktop 
computer. However, it was noticed that even after a 
small number of generations, the correct value of b was 
"winning the race" and so a smaller number of gener- 
ations can be run to determine the interesting regions 
of parameters space for for further exploration. To 
illustrate this. Figure 1101 presents the same result as 
FigureEl using only 200 generations rather than 3500. 
This obviously has a significant computational advan- 
tage. 

Calculating each point in Figures rHand lTITl reauired 
one run of the genetic algorithm, and the population 
was reset each time such that the initial genomes are 
sequences of zeros (corresponding to black pixels). How- 
ever, if each point in these figures are calculated se- 
quentially, the populations need not be reset as the 
solution to the previous point already contains fairly 
good solution, but for a slightly different 6- value. Then 
the GA is given a head start in trying to improve the 
solution, effectively tweaking the previous solution to 
produce a new solution. This results in a dramatic 
smoothing out of the figures (Figure IIH , and allows 
the parameter to be calculated correctly to and accu- 
racy about 3 significant figures in only 200 generations. 

As noted earlier, it should be possible to simply 
generalise this one-dimensional method to include more 
than one free parameter, either by using a large grid 
search (which may take a long time, although it is 
straightforward to devise a parallel computing scheme 
to do this, since each run can be done independently 
of the others), or by using a multidimensional minimi- 
sation method such as Powell's method. 



5 Conclusions 

This paper has introduced a new technique for the 
inversion of gravitational lensed images of extended 
sources. This utilises genetic algorithms to evolve an 
optimal source for a particular gravitational lens model. 
It is seen that this approach successfully recovers the 
source configuration of an idealised gravitational lens 
system. Furthermore, it was demonstrated that this 
genetic algorithmic approach successfully recovers the 
source profile in the presence of noise and can be incor- 
porated into more general gravitational lens optimisa- 
tion schemes. This initial investigation has considered 
only a simple model of gravitational lensing, neglect- 
ing detailed aspects of true gravitational lens systems, 
such as various sources of noise and image smearing 
due to instrumental and atmospheric effects. How- 
ever, due to the forward mapping of this approach, 
these can be added in a straight forward fashion, pro- 
viding an inversion technique that can be applied to 
observed gravitational lens systems. Due to the lim- 
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Figure 10: Best error after 200 generations as a 
function of the value of b. Note that the evolution 
is proceeding at a faster rate, even at this early 
stage, if the lens parameters are correct. 



ited time- frame of this initial project, these aspects of 
the algorithm will be left as further work. 

The most time consuming part of the genetic al- 
gorithmic approach is the calculation of the fitness of 
each member in a generation, scaling with the size of 
the population. For a particular genome, however, the 
calculation of the fitness is independent of the other 
members of the generation. This leads to a simple par- 
allelisation of the approach, with the fitness calcula- 
tion farmed out to individual processors. Furthermore, 
genetic evolution can be driven harder via the inclu- 
sion of parasitic organisms or 'black sheep' ijBobingeJ 
l200Cf) . speeding up the evolution of the genome to fit- 
ter solutions and preventing evolutionary stagnation; 
these too will be incorporated into fuller version of this 
inversion technique. 
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